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Abstract 



Using the tools of noncommutative geometry we calculate the distances 
between the points of a lattice on which the usual discretized Dirac operator 
has been defined. We find that these distances do not have the expected 
behaviour, revealing that from the metric point of view the lattice does not 
look at all as a set of points sitting on the continuum manifold. We thus 
have an additional criterion for the choice of the discretization of the Dirac 
operator. 
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1 Introduction 



Under the generic name of Non-Commutative Geometry JTJ] nowadays is meant 
an approach to the geometry of topological spaces (or their generalization) based 
on the algebra of continuous functions with complex values, or a generalization 
to a non-abelian algebra. The stress thus shifts from the topological space as 
a set of points, to the functions which can be defined on it. While the name of 
noncommutative geometry should in principle only be given to the generalizations, 
it is customary the use of the term even for the case of topological (Hausdorff) 
spaces (see for example [jl], 0, |3| and references therein). In the philosophy of Non- 
Commutative Geometry the main object to study is an algebra A, from which (in 
the abelian case) one can reconstruct an Hausdorff topological space M. It is also 
possible to introduce algebraic tools to analyze metric, differential and algebraic 
topological features. In particular, with the help of a Hilbert Space on which the 
algebra is realized as operator algebra, and of a Dirac Operator, is then possible 
to consider also the metric properties of the space, such as the distances between 
points. 

In this paper we will apply these ideas to the lattice discretizations of man- 
ifolds. The discretization of a manifold, and consequently the change of differ- 
ential equations into finite difference ones, is not always harmless. It is known 
[f|] that not all solutions of the discretized equations in the continuum limit give 
rise to solutions of the related differential equations. In the lattice discretization 
of field theory the well known undesired feature is the appearance of spurious 
states (fermion doubling) (see for example || and references therein). Here we 
will concentrate on the metric properties of the lattice, as derived from the tools 
of noncommutative geometry^ and in particular from the usual discretized Dirac 
operator. 

We find some surprising results. While in the limit of large separation (as 
compared to the lattice spacing) distances behave as expected, for points which 
are only a few lattice sites apart, some interesting anomalies appear. In the case 
of an hypercubic lattice, we find that points belonging to a straight line parallel 
to a coordinate axis ar not 'on a line', in the sense that distances between sites 
two apart is not twice the distance between prime neighbours, and in general 
the distances are not additive, although they become so in the limit of large 
separations. 

tFor a review of approaches to lattice physics based on noncommutative geometry see M and 
references therein. 
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2 Topological and Metric aspects from Algebras. 



The central idea behind non-commutative geometry is in the work of Gel'fand and 
Naimark (see for example |[], which finds a one to one correspondence between 
Hausdorff topological spaces and commutative C*-algebras, that is a commutative 
algebra A on which a norm || • || and an antilinear involution * are defined, such 
that || / || = || f* II, II f*f 11 = 11 /* II II / II and (fg)* = g*f* for f,g G A. The 
algebra A is also assumed to be complete in the given norm. 

The correspondence in one sense is simply made by considering the algebra of 
complex functions on the topological space. In the other direction the correspon- 
dence is between the algebra and the space of its maximal ideals^ A. An ideal / 
is a subalgebra with the property^: ab G / , Va G / ,b G A. A maximal ideal is an 
ideal which is not contained in any other ideal (apart from the whole algebra ^4). 
The space A is equipped with the topology defined by giving the closure of any 
given set of ideals W. The closure is defined as follows 

W = {I £A:lD[)J,JeW} . (1) 

Once the topological space has been reconstructed, one can see that the maxi- 
mal ideals correspond to the functions vanishing in one point, and it is not difficult 
to see that the topology above defined corresponds to the topology of the space 
with respect to which the functions are continuous, 

The power of the method is its ability to be generalizable to noncommutative 
C*-algebras. In this case one cannot in general talk of a structure space, although 
it is still possible in some cases to reconstruct a topological space, which however 
will not be a Hausdorff separable space ||. 

A fundamental ingredient in Connes' noncommutative geometry is a fc-cycle, it 
contains in fact the information about the metric. A k-cycle is a triple (A, H, D) 
consisting of an algebra A, an Hilbert space H on which the elements of A are 
realized as linear operators, and a Dirac operator D which is a selfadjoint operator 
with compact resolvent such that, for a dense subset of A, the commutator [D, f] 
is a bounded operator. For the purpose of this paper, what we are interested in is 
mainly the distance, which can be reconstructed from the /c-cycle in the following 

■^An equivalent reconstruction can be made by considering the space of irreducible represen- 
tations of the Algebra . 

§For noncommutative algebras a distinction between right and left ideals must be made. 
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way: the distance d(x, y) between two points x and y is denned as 



d(x,y)=sup{\f(x)-f(y)\:feA,\\[D,f}\\<l}. , (2) 
/ 

where the norm among operators on H is the one given by the supremum over 
normalized vectors of H. This distance is to be compared with the usual definition 
of (geodesic) distance for a Riemannian manifold with euclidean metric g^ u : 

d(x,y) = inf U,{x,y) (3) 

lry(x, y) being the length of the path 7 from x to y. 

To see the relation between these two definitions of distance, consider the usual 
dirac operator D = ij^^d^ + A^) acting on some dense domain of the Hilbert space 
H of square integrable spinors ^f(x) defined in a IR 4 (A^ is the connection for some 
gauge group G and the spinors belong to some representation of G). Let A be 
the algebra of continuous functions of IR 4 vanishing at infinity, and consider them 
as operators on H acting by pointwise multiplication. The commutator [D,f], 
acting on the spinor is 

[D,/]* = (i^dj^ + if^d^-if^d^ (4) 
= i(Yd»m. (5) 



The norm of the commutator is thus the maximum of yd> 1 f*d fl f, which is also 
equal || to the Lipschitz norm of / defined as 

11 f I, !/(*)- /(y)| fR . 

|| / \\lAp.= SUp r . (6) 

x+y d(x,y) 

Thus 

\f(*)-f(v)\ 



|| [D,f] ||= max Jd»f*dJ =|| / |U ip := sup uv ' (7) 
v x ^ y d{x,y) 

To check that the definitions of distance in (0) and (0) agree, consider that the 
condition on the norm of [D, /] in (g), and equation (|7|) imply 

d(x,y) < d(x,y) . (8) 

Considering now f(t) = d(x, t) as a function of t (suitably regularized), it satisfies 

|| [D, f(t)] \\ = \\d(x,t) \\ Lip =l (9) 

which gives 

d(x,y) = d(x,y). (10) 



4 



Once the notion of distance between pairs of points is given, the other metric 
properties can be obtained with the traditional techniques. 

Thus, we have seen that, in order to reconstruct the topological and metric 
aspects of a topological space one should be given an abelian algebra A realized as 
linear operators on a Hilbert space H, together with another operator, the Dirac 
operator D. 

3 Abelian Discretization of a Manifold 

The discretization of a continuous manifold via a regular lattice is a well known 
technique, which in general provides some simplifications of the problem at hand, 
such as the substitution of differential equations with finite difference equations, 
or the possibility to encode with a finite set of numbers the state of a system. 
Invariably something is missing when a continuous manifold is approximated by 
a lattice. A prime victim is topology, the points of a lattice are equipped with the 
discrete topology, each single point is both open and closed, and every map is a 
continuous function, that is, the inverse image of an open set is open. This for the 
simple reason that every collection of points of a lattice is open. The programme 
of |5|, [TT| is in fact to define generalizations of lattices with non trivial topology. 
In || we have discussed as a lattice with non trivial topology is, in some sense, a 
non-commutative manifold. In this paper we will however limit ourselves to the 
customary, 'abelian', lattices. 

In order to discretize a manifold one chooses a finite or countable set of points, 
so that the state of the system, which in general will be a function on the manifold, 
is now the finite or countable set of values of the function on the points of the 
lattice. Similarly, differential operators are approximated by finite differences 
ones. The hope is that even a limited number of points will capture the essential 
features of the continuous system. Obviously the discretized system cannot have 
exactly the same behaviour of the continuous one. The idea is that by increasing 
the number of points, in other words by taking the continuum limit of the system, 
the correspondence becomes exact. 

From the noncommutative geometry point of view, for a N point discretization, 
the algebra consists of ordered iV-ples of complex numbers, with product and sum 
acting component by component. The number N can be taken to be infinity 
without troubles. From this algebra we reconstruct the abelian lattice, that is a 
set of points with the discrete topology, each point is both closed and open, and 
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therefore so is any subset. On the lattice we then consider the customary discrete 
version of Dirac operator. We then calculate the distances between the points of 
the lattice, with the methods of noncommutative geometry. 

Let us consider in detail the case of a four-dimensional euclidean space-time 
H 4 (with the usual metric). A possible discretization of the manifold IR 4 is then 
accomplished by considering an hypercubic lattice of points with integer coordi- 
nates (in units of the lattice spacing a): 

x = a(ii,i 2 ,i3,U) ,veZ , (11) 

together with the Hilbert space of "square integrable" spinors: 

* x : 5>£**= E*te**«<oo . (12) 

x x,a,i 

where % and a denote spinor and gauge indices respectively. This is the Hilbert 
space commonly called I 2 . On this space the elements of the algebra / = {f x } 
are realized (as operators) as diagonal matrices (we use the same symbols for the 
operator and the element of the algebra): 

fx,i,a;y,j,/3 = fx^xy^a^ij HO SUH1 On X. (13) 

Given this discretization the next question is the choice of the finite difference 
operator which corresponds to the Dirac operator. This choice is not unique: the 
simplest one is that made by Wilson: 

(D^) x = J2r(D^) x , (14) 
where is the covariant discrete derivative in the direction /i: 

In the latter equation x + /i denotes the site displaced from x by one unit along 
direction \i and is the gauge group element (concretely a unitary matrix) 
associated with the link joining x with x + /i, the so-called link- variable. By 
definition, J/"" = (Ug)' 1 . 

Despite its simple look, this Dirac operator hides subtle problems. As is well 
known it is affected by the pathology of "species doubling": in the continuum 
limit, a — > 0, it describes 16 Dirac fermions with degenerate masses. For the time 
being we will not worry about this question and will come back to it at the end 
of section 5. We will now show instead that at the even more fundamental level 
of the definition of distances, an important pathology happens. 
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4 Distances on the Hypercubic lattice 



In this section, we will compute the distance between the points of the lattice 



(|iT|) , using as Dirac operator the finite difference operator (fL4|) . For simplicity, we 
will limit ourselves to the distance dk between the origin of the lattice and a point 
x = (k, 0, 0, 0) of the x— axis (it will be clear from the proof that the distances 
are invariant under translations and rotations that leave the lattice invariant). 
Naively, one would expect to find for dk the following result: 

dk = a\k\ . (16) 

As we shall see in the sequel, this expectation is incorrect. 

The computation of dk consists of two parts: first, we show how to reduce 
the problem to that of finding the distances on a one-dimensional lattice endowed 
with a "Dirac" operator given by the flat discrete derivative in one dimension, 
namely eqn. (|15|) with fi = 1 and link-variables all equal to I: 

pi 0) ^ = ^(^, + i-^,-i) • (IT) 

We then proceed in the next section with the computation of the distances for 
this simpler one-dimensional lattice. 

In order to achieve the reduction to a one- dimensional problem, we start by 
proving that the supremum of (2) can be attained considering functions which 
depend only on the x— coordinate of the lattice sites. In fact, let / = f(ii,i2, iz^ h) 
be any function saturating eq.(2): 

II [DJ] \\<l, 

d k = \f(K 0,0,0) - /(0,0,0,0)| . (18) 
Consider now the function f( x ' defined as: 

/^(ii.*2,i3,*4) = /(ii, 0,0,0). (19) 

That is we consider a function which coincides with / on the X\ axis, and is 
constant along the other three axis. By construction f( x > depends only on i 1 and 
gives for dk the same result as /. Then, if we show that || [D,f( x >] ||< 1, the 
statement will be proven. Now, for any normalized \1/ G I 2 , we have: 

II [DJ^ || 2 =|| 7 1 ^i,/ (x) ]^ f=ll [Di, f (x) }* f • (20) 



7 



The norm || [Di, fw] || is easily seen to be independent on the link variables. 
Consider in fact the following gauge transformation U: 

if x = {m) x = g x ^ x (21) 

where g x are unitary matrices such that: 

UxQx+i = U~\g x -i for i x > 

UZ-i9x-i = U l x g x+1 for ix < 1 . (22) 
It is easy to verify that: 

|| [D u f&] ||= sup || [£>!,/<*>]* ||= sup || [D u fM\m || = 
H*ll 2 =i ll*ll 2 =i 

sup || [£>?\ /<»>]¥ || = || [Df\fW] || , (23) 
ll*ll 2 =i 

Now, the Hilbert space I 2 can be decomposed into a direct sum of subspaces left 
invariant by the operator [D±*\ f^\. 

I 2 = #(<>W) , (24) 

l,m,n;j 

where 

H (i,m,n-,j) = ^ . = g un i ess ^ _ ^ i 3 = m , i A = n and i = j} . (25) 

Moreover, the restrictions of [d[ , f^] to these subspaces are all equal to each 
other and then we can write: 

n rn m(0) i-fxli i, 

\(l,m,n;j) 



[D?\fM] ||= sup {|| [D?>,fW] |k m , n;j) } 



(l,m,n;j) 

= 11 [M 0) ,/ W ] ||(0,0,0;l) =11 [DxjM] ||(0A0;l) , (26) 

where || [Di,/^] ||(j, m ,ny) denotes the norm of the restriction of [Di, f^°>] to 

But it is now easy to check that, for any normalized \1/ e f/wWi^ we h ave: 

11^1,/^ll^lltA/]*!^!. (27) 
and then, in view of (^) and (|26|), it follows that: 

ll[£>i,/ W ]IHI[^ 0) ,/ W ]ll(o.o^i)<l. (28) 

What this equation proves is that the distances dk coincide with those for a one- 
dimensional lattice endowed with the flat-derivative operator of eqn. re- 
garded as Dirac operator. We have thus reduced the calculation of the distances 
on a one-dimensional lattice. We calculate these distances in the next section. 
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5 Distances on the one-dimensional lattice 



In this section, we compute the distances on a one- dimensional lattice: 



Xk = ak ,k 62 . 



(29) 



with the Dirac operator given by eqn. (fT^). For economy of notation this oper- 
ator will be denoted in this section by the letter D. First of all, it is clear that 
distances are translationally invariant. It is enough then to compute the distance 
dt between sites xq and Xk- Our proof goes as follows: in order to obtain the 
supremum in (Q) we first show that this can be attained considering only real 
'step functions', functions which are constant except for a finite set of points, on 
which they monotonically increase. The distance between the points is then sim- 
ply the total step. With these functions the problem becomes the one of finding 
the maximum step the function can have, subject to the restriction that the norm 
of the commutator is less than one. For step functions this last requirement will 
mean that the eigenvalues of a finite matrix are all less than one in modulus. 

As a first step, we will prove the distance dk can be obtained by taking the 
supremum of (2) over the set of real functions = {// ; i G Z} such that: 



In fact, let / = {fj ; j G 2} be an arbitrary complex function on the lattice. We 
then have: 



In what follows we will set 2a = 1, the changes for a general a being trivial. For 
any \1/ G h, consider a new vector \1/ such that: 



Such a vector always exists (for example, one can choose \I/o = and = 
and use eq. ([H]) to recursively determine ijjj for the other j's). We then have: 




(30) 



arg((/ i+ i - fj)^ j+ i) = arg((£ - Vj G TL. 



(31) 
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< E(l/i+i-/iH^ + il + l/i-/i-ill*i-il) 2 =H [A/]*H 2 (32) 



j 
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It follows that 



|| [D, f] || 2 = sup - /,-|c i+1 + |/, - / i+a |c,_ 1 ) a . (33) 

Cj eIR + ,^. c 2=i j 

Consider now the real, monotonic function F defined by: 

F J+1 -F 3 = \f 3+1 -f 3 \ , F = 0. (34) 
Considerations similar to those applied to / imply that the norm of [D, F] is again 



given by equation (|33|). It follows then that 

|| [DJ] || = || [D,F] || . (35) 

Moreover, we clearly have 

fc-i fe-i 

\fk — fo\ = I ^2 Um+l — fm)\ < ^2 \fm+l — /m| = F k — Fq. (36) 
m=0 m=0 

Finally, consider the function which coincides with F for < i < k. Use of 
equation (^) implies that 

|| [D,fM] \\<\\ [D,F] || . (37) 

From the last inequalities it then follows that in order to compute the distance 
between any two points we can restrict the supremum to the functions 

4 = sup{/f - fi k) : || [£>,/«] || <1} . (38) 
/W 

Making use of (|T7|), it is now easy to check that the norm of [D, f^] is equal 
to the maximum eigenvalue r of the square, symmetric and real (k + 1) x (k + 1) 
matrix H: 

H m ,n = ^■m$m+l,n + ^m-l^m-l,ra TTI = 2, ■ ■ ■ , k] Tl = 1, • • • , k + 1. 

Hi >n = Ai5 n)2 ; Hk+i, n = Afe5 n) fc , (39) 

where A/ = — /^{, i = 1, • • • , k. Then equation (|38|) for the distance can be 
rewritten as: 

f fc "I 

(40) 



4 = sup I V Aj : r < 1 L 

/« U=i J 



Now, according to a theorem of linear algebra on matrices with non negative 
elements |TT[, the maximum eigenvalue r of H is less than 1 if and only if all the 
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leading principal minors if n (A 1; • ■ ■ , A n _i), n — 1, • • • , k + 1 of the matrix I — 
are positive: 



Fx = 1 > ; H 2 = 1 - A\ > ; • • • = det(K - H) > 0. (41) 

We notice the following recursive relation among the H n : 

H n+1 = H n - A=ff n _! , 1< n < fc + 1 . (42) 

The determination of the distances is thus reduced to finding the supremum of 
E*=i Ai in the open region X of IR fc specified by the set ([|l]) of inequalities. It is 
clear now that J2i=i Aj cannot have a local maximum in the interior of X. It then 
follows that the supremum must be reached at some point of the frontier dX of 
X, which is the set of points y such that: 



1. At least one among H n (y) is equal to 0. 

2. There is a sequence of points of X converging to y. 

Now, dX = Uf=2 9X h where dXi is the open surface given by the equation^ 

^^{(A^-.-.Afc) : A m >0; H j > 0, j = 1, • • -i - 1; H, = 0}. (43) 

Then dk = max{Aj, i = 2, • ■ • , k + 1}, where Aj is the maximum, when it 
exists, of SjLiAj on dXi. We can now prove that, for i < fc + l, it must be 



Aj = <ij_i+c4-j- In fact, let (A^, ■ • • , A^') be the point of dXi such that max(Ai + 
■+A k ) = (A[ i] + - ■ ■ + A i k i) ) = Aj. Then, it follows from the the recursion relation 



.(0 



that Af> = 0. In fact, if Af } ^ 0, we would have H i+1 = Hi- Af'H^ < and 



then (Ai \ ■ ■ ■ , A^ iJ ) ^ <9X. But Af ' = implies that the matrix H block diagonal 
and this obviously implies the statement. Then, if the distances di, i — 1, • • • , k— 1 
are known we need only compute \k+i- The latter can be found with the method 
of the Lagrange multipliers: one looks for the solutions of the following system of 
algebraic equations in the unknowns (A^ +1 \ • • • , Aj^ +1 ^) and fi: 



.0 







OA, 



£ A + 



i=l 





fe+i = o. 
(fc+lh 



(44) 



(a posteriori one has to check that H^A^ , • • ■ , A^ 1,1 ) > 0, z = 2, • • ■ , fc.). The 
first few distances can easily be computed analytically. For an arbitrary a one 

^On dXi, the recursion relations imply Hi = 0, for I > i, if i/,; = 0. 
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Sites Apart 


1 2 3 


4 


5 


6 7 


8 9 


10 


11 


12 


13 


Distance 


2 2.83 4 


4.90 


6 


6.93 8 


8.94 10 


10.95 


12 


12.96 


14 



Table 1: The distances calculated analytically for the site differences less than 4, 
numerically otherwise. 



gets: 

di = 2a 

d 2 = 2\[2a 

c?3 = 4a 

d 4 = 2y/Ea. (45) 

The distances dk, for k > 4, have been determined numerically. In the table the 
distances for sites up to 13 apart are shown. The results show that the sites with 
k even behave differently from those with k odd. The result can be expressed by 
the following formulae: 

d 2m = 2(/ + l)a 
d 2l = (2/+l)a-e 2/ (46) 

where e 2 \ > 0. We thus see that there is an overall constant "anomaly" equal to a 
for all sites, plus an extra one, e 2 i, for the even sites, which vanishes for 21 — ► oo. 
The conclusion of this analysis is that, form the point of view of the finite difference 
operator ([14]), the lattice does not look at all as a subset of equidistant points of 
the real line. Preliminary results seem to indicate that in order to get the desired 
distances, one is forced to use non-local difference operators. 

As pointed out at the end of section 3, the 'naive" discrete Dirac operator 
( |T4D is plagued with the problem of species doubling: in the continuum limit it 
describes 16 Dirac fermions with equal masses. Many possible remedies to this 
difficulty have been considered in the literature and one may wonder whether 
they will cure the pathology in the distances as well. We shall see immediately 
that this is not the case. It will be sufficient to consider, as an example, the 
modification of the discrete Dirac operator introduced by Wilson: one adds to 
( |i"4"|) a term proportional to the discretized covariant laplacian, whose effect is to 
make the masses of 15 out of the 16 fermions of the order of the inverse of the 
lattice spacing. In this way the spurious fermions disappear in the continuum 
limit. The modified operator Dy/ reads : 

{Dw*) x = {D*) x + L. E^vD^ + U~^ x ^ - 2* x ) . (47) 
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where r is an arbitrary non vanishing real constant. If we take this as the Dirac 
operator to be used in the computation of the distances instead of flU]) , it is easy 
to verify that the only change in the distances will be an overall rescaling by a 
factor (1 + r 2 ) -1 / 2 (this shows that from the distances point of view the inclusion 
of the Laplacian actually worsens the situation instead of improving it, because 
now we do not recover the correct distances even for large separations). The fact 
is that, as is clearly shown by the computations of section 3 and 4, the distances 
mainly depend on the choice of the discrete flat derivative in one dimension. 

6 Conclusions 

The first conclusion we can draw is that the Dirac operator usually used on the 
lattice does not reproduce the distances one would expect. This opens another 
problem. We have three concepts of distance at hand: the naive one of equation 
([16]); the one we calculated in sections 4 and 5, which however in the limit of large 
separations goes to the naive one; finally we can envisage some sort of 'geodesic' 
distance: 

i 

d(x k ,xi)= J2 d(zi,2C f _i) . (48) 

i=k+l 

This last distance in our case is twice the naive distance for all separations inde- 
pendently of the lattice spacing a, and thus has a different continuum limit. 

One would of course like that all those distances be the same, as in the con- 
tinuum case. Since the only ingredient of this calculation is the choice of the 
discretization of Dirac operator, we think that we are seeing another pathology of 
the choice (fH]). We do not know of a discretization of the Dirac operator which 
reproduces correctly the distances. 

In conclusion, the results of this paper show that the concepts of noncom- 
mutative geometry, apart from representing a fashinating generalization of known 
geometrical notions, may also suggest new points of view about conventional prob- 
lems, like the discretization of a field theory on a lattice. We have seen that if, 
on a lattice, one has assigned distances (inherited, say, from the continuum), non- 
commutative geometry provides additional, stringent criteria for the choice of a 
discretization of the Dirac operator. 
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